Kerr nonlinear coupler and entanglement 
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We discuss a model comprising two coupled nonlinear oscillators (Kerr-like nonlinear coupler) 
with one of them pumped by an external coherent excitation. Applying the method of nonlinear 
quantum scissors we show that the quantum evolution of the coupler can be closed within a finite 
set of n-photon Fock states. Moreover, we show that the system is able to generate Bell-like states 
and, as a consequence, the coupler discussed behaves as a two-qubit system. We also analyze the 
effects of dissipation on entanglement of formation parametrized by concurrence. 
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I. INTRODUCTION 

Quantum entanglement seems to be one of the most 
striking phenomena of quantum physics. It is not only 
one of the most fundamental concepts of quantum in- 
formation theory, but also allows investigation of many 
features of nonlocal properties of quantum systems 
Various aspects of the entanglement and its generation 
have been discussed in numerous papers, especially from 
the point of view of quantum information applications in- 
cluding quantum key distribution j^i], superdense coding 
0, quantum teleportation fast quantum computa- 
tions entanglement-assisted communication '7'| or 
broadcasting of entanglement Q . 

In this paper we shall concentrate on the dynamics 
of the Kerr nonlinear coupler and its ability to produce 
quantum entangled states. Since the pioneering works of 
Jensen 9| and Maier llOll . the nonlinear couplers attract 
an increasing interest Illll-f T?! (for reviews see [T^ IT^). 
As shown in classical ^9lll0|| and quantum [lli | models, the 
Kerr couplers can exhibit self-trapping, self-modulation 
and self-switching of the energy of the coupled modes. 
These phenomena have potential applications in opti- 
cal communications as, e.g., intensity-dependent routing 
switches. Among various other quantum statistical prop- 
erties, it has been shown that the Kerr couplers can be 
a source of sub-Poissonian and squeezed light *T^-p^. 
Another group of papers concerns the correspondences 
between the quantum and classical dynamics of such sys- 
tems Tlj and their chaotic dynamics, including synchro- 
nization effects fr^ . 

Quantum optical systems based on Kerr nonlinearity 
have been applied for various quantum information pur- 
poses including entanglement purification (2Qj . complete 
quantum teleportation |^| , or realization of qubit phase 
gates l2.2j. Here, we present another simple quantum 
information application Kerr nonlinearities, namely for 
generation of entangled optical qubits from classical light. 

We are interested here in a simple model comprising 
two quantum nonlinear oscillators located inside one cav- 
ity. These oscillators are linearly coupled to each other, 
while one of the oscillators is excited by an external co- 
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FIG. 1: Scheme of a pumped nonlinear coupler, described by 
Hamiltonian Q, implemented by two ring cavities (a and 6) 
filled with Kerr media, where cavity a is being excited by a 
single-mode external classical field. 



herent field of a constant amplitude. For this model we 
shall answer the questions, whether it is possible to close 
the dynamics of the excited nonlinear coupler within a 
finite set of n-photon states and, which is the main sub- 
ject of this paper, whether nonlinear excited coupler can 
be a source of maximally entangled (ME) states. 



II. THE MODEL AND SOLUTIONS 

The model of the Kerr nonlinear coupler discussed here 
contains two nonlinear oscillators linearly coupled to each 
other and, additionally, one of them is coupled to an ex- 
ternal coherent field as presented in figure 1. We as- 
sume, that this excitation is linear and has a constant 
amplitude. This system can be described by the follow- 
ing Hamiltonian: 



H — Hnl + Hint + Hext 



(1) 



where 
Hnl 

Hint 
Hext 



m + H^^l = ^{a^fa^ + f{b^)H', (2) 



ea)b + €*ah\ 
aa) -\- a* a. 
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We see that Hnl describes nonlinear oscillators, Hint 
corresponds to an internal coupling, whereas the term 
Hf,xt describes a linear coupling between the external 
field and the mode of the field inside our cavity corre- 
sponding to the oscillator a. The parameters Xa and 
Xb are nonlinearity constants of the oscillators a and b, 
respectively, e describes the strength of the oscillator - 
oscillator coupling, whereas a is a strength of the exter- 
nal excitation of the oscillator a. It is worth noting that 
our Hamiltonian Hint does not include nonlinear coupling 
proportional to b'^ba^ d but only the linear one, described 
by ©. Nevertheless, the same Hamiltonian Hnl + Hint 
as ours was used, e.g., by Bernstein and Chefles and 
Barnett to describe the nonlinear coupler. 

In the first part of our analysis we neglect damping 
processes in our model, thus the system evolution can 
be described by a time-dependent wave- function. This 
function can be written in the n-photon Fock basis as: 



dt 



(7) 



m)) 



(t)|n)a|m)b 



(5) 



n,m— 



where Cn,m{t) is a complex probability amplitude of find- 
ing our system in the n-photon and m-photon states for 
the mode a and b, respectively. 

We have included here an external coupling and there- 
fore, the energy inside the cavity is not conserved. As a 
consequence, we can expect that in the evolution of the 
system many of the states corresponding to great number 
of photons will be involved. However, we can overcome 
this difficulty by applying the nonlinear quantum scis- 
sors method discussed in \24\ (for the discussion concern- 
ing quantum states defined in finite-dimensional Hilbert 
spaces and the methods of their generation see the re- 
view papers and the references cited therein). 
Namely, it is seen from the form of Hn l that this Hamil- 
tonian produces degenerate levels of the energy equal to 
zero, corresponding to the following four states: |0)a|0)b, 
|l)a|0)6, |0)a|l)b and |l)a|l)6- Morcover, all couplings 
discussed here have constant envelopes, and similarly as 
in (24 |. we assume that they are weak. Therefore, we can 
treat transitions within the mentioned set of the states 
as that of resonant nature. The evolution of discussed 
system is closed within the set of these four states and 
interactions with other states can be neglected in our 
approximation. Thus, the wave-function describing our 
model can be written in the form: 

IV'(t)) = Co^o(i)|0)a|0)6 + Ci,o(<)|l)a|0)h 

+ Co,l{t)\0)a\l)b + ClAt)\l)a\l)b (6) 

and hence, the equations of motion for the system are: 



dt 
. d 

dt 
■ d 

*:i:Co,i = e ci,o + a ci 1, 
dt 



To solve these equations we need to find roots of the 
forth-ordcr polynomial, which leads to a very compli- 
cated and unreadable form of final formulas. Therefore, 
although it is formally possible, we shall not write the 
analytical solutions in their most general form and we 
will restrict our considerations to the case of real a — e. 
Moreover, we assume that for the time t — the both 
oscillators were in vacuum states, i.e., 



mt = 0)) = \0)a\0)b. 



(8) 



Then we get the following solutions for the probability 
amplitudes Q.j (i, j = 0, 1): 

coo(^) = cos(a;t) cos(?/i) H — ^ sm{xt) sm{yt) , 
2 

ci.o(i) = —i—^cos{xt)sm{yt), 



Co At) = -—^sin{xt)sm{yt), 

ci i{t) = i[^= cos(xt) sin(yt) — sin(xt) cos(2/i)] (9) 

where x — a/2 and y — \/hx. The solution ^ is valid 
under the condition e = a (j = a, 6), which implies 

that it is apparently independent of nonlinearities Xj- 
But it should be stressed that the corresponding nonlin- 
ear Hamiltonian H^Li given by |2Jl, is responsible for the 
truncation of the infinite-dimensional state to the finite 
superposition, given by ©. Otherwise, if Xj were not 
much stronger than e and a, the state generated would 
not be truncated to the finite superposition (jSJ and the 
probability amplitudes Cn,m(t) would depend explicitly 
on nonlinearities Xj- 

To check our solutions we can calculate the probabil- 
ity amplitudes numerically in a basis expanded to the 
states corresponding to greater number of photons than 
discussed here (for the model discussed our considera- 
tions are restricted by the resonances to the one-photon 
and vacuum states only). We perform the calculations 
following the method discussed in 26], and first we con- 
struct the unitary evolution operator U applying the full 
Hamiltonian shown in (Q: 



U 



exp 



(-iHt^ 



(10) 



Then we are able to obtain the wave-function li^it)) by 
acting the operator U on the initial state of the system, 
and for the case discussed here we have: 



- Uit)\0)a\0)b. 



(11) 



Figure 2 shows both analytical and numerical results of 
our calculations. We see very good agreement between 
these two methods, so the model based on the resonances 
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FIG. 2: Probabilities for finding tfie coupler in the |0)a|0)6 
(dotted curve), |l)a|0)6 (dashed-dotted curve) |0)a| 1)6 (dashed 
curve) and |l)a|l)i, (solid curve) states from the analytical re- 
sults and their numerical counterparts (cross marks). The 
nonlinearity coefficients Xa ~ Xb = 25 and the coupling 
strengths e = q = 7r/25. 




FIG. 3: The same as in figure 2 but for the states: (a) |0)a|2)6 
and (b) \l)a\2)b. 

works very well. Moreover, from the numerical results 
for the probabilities corresponding to the states |0)a|2)b 
(figure 3a) and |l)a|2)fc (figure 3b), we see that the states 
corresponding to the numbers of photons higher than one 
are practically unpopulated. It is worth mentioning that 
our numerical calculations have been performed in the m- 
dimensional Fock basis, where to ~ 20 for each subspace 
associated with a single mode of the field. 

III. COUPLER AND ENTANGLEMENT 

The time-evolution of the probability amplitudes can 
give some information concerning entanglement in our 
system too. For instance, if we see in a figure that prob- 
ability corresponding to one of the discussed states is 
equal to another one and additionally, the both are equal 
to 1/2 (|cijp = |cfc,iP — 0.5 for every i,j,k,l), we know 
that the system generates ME states. Obviously, this 
method of finding entangled states is not very accurate, 
especially for the case when we should observe and com- 
pare various and often rapidly oscillating probabilities. 
Therefore, we apply another method convenient for find- 



FIG. 4: Probabilities for finding the coupler in the Bell- 
like states. The nonlinearity coefficients Xa and Xb, and the 
coupling strengths e and a are identical to those of figures 
2-3. 

ing entanglement in the system. Namely we shall express 
obtained wave-function in the Bell basis: 

\^) = h\Bi) + 62IB2) + 53IB3) + h\B4) (12) 

where the states \Bi), i = 1,2,3,4 are Bell-like states 
that can be expressed as functions of discussed here n- 
photon states (Bell-like states differ from the commonly 
discussed Bell states in the existence of the phase factor 
- for the case discussed here, one of the n-photon states 
is multiplied by z) : 



\Bi) 


- ^(|l>a|l>. 


+ *|0)a|0)fc). 






+ i\l)a\l)b), 


\B3) 




-*|l>a|0)fc). 


\B4) 


= ^(|l>a|0)6 


-*|0).|l)fc). 



(13) 



These states are maximally entangled states, and there- 
fore, for the cases when = 1, i = 1,2,3,4 also our 
system evolves into ME state. Figure 4 shows probabili- 
ties corresponding to the Bell-like states as a function of 
time. Moreover, all parameters describing our system are 
identical to those of figure 2. We see, that for the time 
t ~ 115 we get the state \Bi) and for i ~ 80 the state 
|i?2) is generated with high accuracy - our system be- 
comes maximally entangled. This entanglement involves 
the states |0)a|Of,) and |l)Q|lb). Of course, one should 
keep in mind that plots in figure 4 are for the proba- 
bilities not for their complex amplitudes and hence, we 
get the Bell-like states from (|13|l with some phase fac- 
tor. Nevertheless, our states are maximally entangled. 
Moreover, figure 4 shows that the values of probabilities 
for the states {B^} and \B4) can reach maximally 0.8. 
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As a consequence, the states |l)a|Oh) and |0)a|l6) can- 
not be maximally entangled for the initial vacuum states 
|0)a|Ob). But generation of IB3) and I-B4) would be possi- 
ble by assuming that the system is initially in the states 

|l)a|Of,) or \0)a\h). 

The Bell-like states are maximally entangles, how- 
ever, they are not the only entangled states that could 
be produced by the system. Therefore, to measure the 
entanglement degree of the system we apply the mea- 
sure that is referred to as the concurrence. This quantity 
proposed by Wootters |23| is one of the most commonly 
applied measures of the entanglement. The concurrence 
for two qubit states is defined as 

C = max {0, Ai - A2 - A3 - A4} (14) 

where Ai (z = 1, . . . ,4) are the square roots of the eigen- 
values of the matrix 

p = p {a'; al) p* {a'; ® a^) (15) 

where cTy"'''^ are Pauli matrices defined in subspaces cor- 
responding to the modes {a, 6}, and the eigenvalues Ai 
appearing in (|14|l should be taken in decreasing order. 
Concurrence takes values from to 1, where for unen- 
tangled states it vanishes, whereas for ME states it is 
equal to one. 

Damping is the main and unavoidable source of deco- 
herence which can easily destroy entangled states. Hence, 
for our results to be applicable in real physical systems, 
we present a numerical analysis of the damping effects on 
the concurrence. Let us assume that the leakage of pho- 
tons from the cavities a and b is described by the rates 
Ka and Kb, respectively. Starting from Hamiltonian 
and defining the collapse operators by Cq — -v/Sk^o and 
Cb = ^/2Kbb, we can write the time-independent Liouvil- 
lian in the standard Linblad form 




(16) 

The evolution of the density matrix p{t) in the dissipative 
system can be found numerically as a series of complex 
exponentials exp((Tfct) given in terms of the eigenvalues 
CTfe of the Liouvillian £, given by (|16|l . 

Figure 5 shows the plot of the concurrence evolution 
for the case discussed here - we have single external ex- 
citation of the coupler and we assume that all couplings 
existing in the system are weak {a — e — tt/25) in com- 
parison to Kerr nonlinearities Xa — Xb = 25. Various 
curves in figure 5 correspond to concurrence evolutions 
with different dissipation rates. We see that the varying 
in time concurrence is modulated by an oscillation of low 
frequency. As a consequence, several maxima appearing 
here are of various values. Two of them for dissipation- 
free evolution (depicted by solid curve), which are the 
closest to unity, correspond to the formation of Bell-like 
states \Bi) and IB2) discussed earlier and shown in fig- 
ure 4. As a consequence, we can treat our system as a 




Time 

FIG. 5: Concurrence for the excited nonlinear coupler for 
various cavity leakage rates Ka = Kb equal to (solid curve), 
10"" (dashed curve) and 10 ^ (dot-dashed curve). The quan- 
tities Xa, Xb, e and a are identical to those of figure 4. 



source of ME states for low dissipation. On the scale 
of figure 5 and for the chosen coupling parameters, the 
differences between the evolution with the leakage rates 
Ka = Kb < 10~^ and the dissipation-free evolution are 
invisible. However, higher leakage rates beyond short 
time evolution cause essential deterioration of the con- 
currence limiting the effective generation of ME states. 
Thus, the results of our calculations indicate that the sys- 
tem discussed is highly sensitive to the dissipation pro- 
cesses. Even relatively small losses from the cavity are 
able to destroy the entanglement. Therefore, we should 
assume that we deal here with a very high Q cavity that 
is capable to preserve practically the whole radiation field 
located inside. However, this assumption is very desir- 
able from our point of view. For this case the coupler 
can be weakly excited by external fields only - the less 
photons can escape from the cavity through the mirror, 
the smaller number of photons can be injected inside this 
way. 



IV. CONCLUSIONS 

In this paper we have discussed a model of nonlinear 
coupler linearly excited by a single-mode coherent field. 
We have shown that the evolution of the system is closed 
within a finite set of states and only |i)a|i)b = 0, 1) 
states are populated. We have applied here the method 
used for the nonlinear quantum scissors and have 
found some analytical formulas for the probability ampli- 
tudes corresponding to these states. We have shown that 
starting from the vacuum state |0)a|0)f, of our system, its 
evolution leads to Bell-like states generation. Moreover, 
we have calculated the concurrence and its behaviour in- 
dicates that the ME states are produced for the system 
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if the photon leakage rates out of the cavities are less 
than 10"^ for the chosen coupling and nonlinearity pa- 
rameters. Moreover, we have shown that the concurrence 
exhibits some modulation effect as a result of the exis- 
tence of various couplings in our system. For each of the 
couplings we have some frequency and their interference 
leads to some long-frequency oscillations in the system. 

We see, that our model, despite its simplicity, exhibits 
intriguing features. We can say that the properties of the 
system discussed here are much desired from the point of 
view of the physical properties of the nonlinear couplers. 
Our scheme can be used for generation of the entangled 
optical qubits from classical light, which is a basic but 
rather simple quantum information problem. Introduc- 
tion of conditional measurement in this scheme, along 
the lines of Duan et al. proposal 20] , is probably worth 
of further study from the point of view of more sophisti- 
cated quantum information applications |20l l2lL . 

Finally, we mention about the experimental feasibil- 
ity of the presented scheme. Since our solution is ap- 
plicable only when the cavity-field intensities are very 



small, so an objection arises that the Kerr nonlineari- 
ties are usually negligible in this case. However, the re- 
cent breakthrough advances on nonlinear optics involv- 
ing very weak light fields show that the nonlinearities can 
be enhanced by several orders of magnitude in ultracold 
atomic systems using electromagnetically induced trans- 
parency when resonant optical absorption is eliminated 
f j28j and references therein). In particular, giant Kerr 
nonlinearities have been theoretically predicted (2^ and 
first experimentally measured to be ^ 10^ greater than 
those in the conventional optical materials [s^l- Thus, 
we believe that the scheme discussed here can be feasible 
experimentally. 
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